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P  pressure,  psfa 

P.  stagnation  pressure,  psfa 

?"'  pressure  at  Mach  Number  unity,  psfa 

PB  back  pressure,  psfa 

P#  exit  pressure,  psfa 

P_  pressure  before  the  shock,  psfa 
x 

P-  pressure  behind  the  shock,  psfa 

PR  pressure  ratio:   P2A1 

<4  constant  heat  flux,  Btu/lbm-ft 

R  gas  constant,  f  t-lbf/lbm- *R 

S  entropy,  Btu/lbm-  *R 

T  temperature,  °R 

T#  temperature  at  Hach  Number  unity,  *R 

T.  stagnation  temperature,  °R 

Tn  mean  temperature,    *R 

Tr  temperature  ratiot  ^^An 

Tg  temperature  where  the  pressure  is  the  back  pressure,  CR 

Tx  temperature  before  the  shock,  °R 

Ty  temperature  behind  the  shock,  *R 

V  gas  velocity,  fps 

V*  gas  velocity  at  Mach  Number  unity,  fps 

VR  gas  velocity  ratio,  Vg/V^ 

Vx  gas  velocity  before  the  shock,  fps 

Vy  gas  velocity  behind  the  shock,  fps 

W  weight  of  gas,  lbj> 


V*  the  weight  of  gas  between  any  section  to  the  section  at 

Mach  Kuaber  unity,  Ibf 

2  elevation,  ft 

!•  the  elevation  where  Maoh  Number  la  unity,  ft 
Zrfmx   the  maxiun  elevation,  ft 

°<  1,  °^2,  °^3#  °fy»    parameters  defined  on  page  17 

Pi.  |32,  (33t  (BUf 

P5,  P6,  @7,  ^8,    parameters  defined  on  page  31 

^1,  ^2,  <*3,  ^k,    parameters  defined  on  page  92 

0  density,  lbg/ft3 

p "'  density  at  Maoh  Kuaber  unity,  lbg/ft"* 

fx  density  before  the  shock,  lb^/ft^ 

pj  density  behind  the  shook,  ibg/ft^ 

?U   density  ratio:   Pg/pi 

2 

f    frictlonal  stress,  lb^/ft 


INTRODUCTION 

Moat  of  the  problems  encountered  In  gas  dynamics  deal  with 
horizontal  flow.  In  this  report,  vertical  flows  of  a  perfect 
gas  are  considered,  in  which  the  change  in  elevation  is  taken  as 
the  independent  variable.  Three  cases  are  considered:   (1)  re- 
versible adiabatlc  flow,  (?.)   reversible  diabatic  flow  and  (3) 
irreversible  adiabatic  flow.   Particularly  investigated  is  the 
weight  of  gas  existing  within  the  pipe  between  any  two  levels  of 
elevations,  the  conditions  that  produoe  normal  shocks  and  the 
locations  of  these  normal  shocks. 

The  four  fundamental  principles  governing  the  motion  of 
compressible  fluids  are  the  law  of  the  conservation  of  mass, 
Newton's  laws  of  motion,  and  the  first  and  the  second  laws  of 
thermodynamics.  The  restrictions  and  hypothesis  for  the  present 
investigation  are  as  follows: 

(1)  One  dimensional,  compressible,  fluid  flow 

(2)  Steady  flow 

(3)  Perfect  gas 

(ij.)  Vertical  circular  pipe  with  constant  cross-sectional 
area 

(5)  Negligible  boundary  layer  effects 

(6)  Negligible  thickness  of  the  normal  shock. 

Based  on  the  above  principles  and  assumptions,  the  general 
analytic  expressions  are  derived  and  presented  in  the  following 
pages.  Their  results  are  also  illustrated  with  numerical  exam- 
ples. 
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Fig.  1.  Control  surface  for  analysis 
of  isentropic,  constant- area,  vertical  flow, 
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Fig.  2.  Control  surface  for 
analysis  of  reversible, 
diabatic,  constant- area, 
vertical  flow. 
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Fig.   3.   Control   surface  for 
analysis   of  irreversible, 
adiabatic,    constant-area, 
vertical! flow. 


ISENTROPIC  PLOW 
Physical  Equations 
Perfect  Gas  Equation: 

JL.     .  JLL  +  -4L-  (l) 

I      y      T 


Definition  of  Mach  Number: 

dM2    B   dV2        dT 
M2         V2    *   T 


Energy  Equation! 


Equation  of  Continuity: 

dg       1   dV2 
?    +   2   V2 


Momentum  Equation: 

JR.   *    gdZ       VdV 
f         It         So 


Equation  of  Sound  Speed: 


(2) 


CpdT  ♦  JHJ-   ♦  -g|-   •  o  (3) 


<fc> 


(5) 


dC  m       1   ■  ,M 

A    ■    2    I  *  ' 


Equation  of  Imj  ulse  Function: 

dF       dP         kdM2  (7\ 

F        P        1  +  W«2 


The  Change  of  the  Properties 

The  differential  relations  of  the  propertiee  are  established  as 
follows. 

(a)  Pressure: 

Combining  Eq.  (k)   and  Eq.  (5) 

-£|L.  .  -£L.  (  M2  •  i  )  (8)1 

(b)  Temperatures 

Eqs.  (3)  and  Eq.  (5) 

(9) 


dT  - 

dP 

(c) 

Sound  Speed 
Eq.  (6) 

dC   . 

c 

i: 

1  dT 

2  I 

(d) 

Gas  Speed: 

Combining  Eqs,  (1), 

,  (3) 

and 

(8) 

yields 

-f '•{*'•* 

k-l/*,? 
-rr-(Mc« 

•") 

♦  t 

(> 

►  (k-1) 

(•) 

Density: 
Eq.  (I*) 
d? 

1  dV2 

2  V2 

•  0 

(6) 


1 
Developed  in  the  Appendix,  p.  70. 

2 

Developed  in  the  Appendix,  p.  70. 


(f) 

Mach  Number: 

Combining  Eqs, 

(9)  and  (1) 

\£       »   (k 

"  M    ? 

Combining  Eqs. 

(2)  and  (l|.) 

d^2 

,r2 

2  —    - 

dT 
T 

Developed  in  the  Appendix,  p,  71. 

2 

Developed  in  the  Appendix,  p,  72. 


(11) 


(12) 


Substituting  Eqs.  (12)  into  Eqs.  (11) 

J^L     .  .  (  k  +  1  )  4Ji  (12-a) 


(g)  Stagnation  Temperature: 

From  the  definition  of  stagnation  temperature 

T0   «  T  (1  +  -*^  M2) 

dT0     dT       (  k  -  i  )M2   dY2 
T6      T       2  ♦  (k-l)M2    M2 

Substituting  Eq,   (12-a)    into  the  above 

jgs.,    |l  LaJWLU  «2)   -Al  ,,Mi 

TT        2  +  (  k  -  1   )  M2  f  (13) 

(h)  Stagnation  Pressure: 

_d£<L     ,j£     2k(l-M2) (  ,  )2 

P0        f    2  4  (k  -  1)  M2  l**1 


(i)  Impulse  Function: 
dP 


d?   k  (1  -  M2) 
f   "TTkH2 


(15) 


Prom  Eqs.  <fc) ,  (6),  (8),  (9),  (10),  (12),  (13),  <UO  and  (15), 
the  changes  in  the  fluid  properties  are  summarised  for  upward 
and  downward  flows,  and  for  subsonic  and  supersonic  velocities. 


dZ  ?0 
(Upward  Plow) 

dZ<0 
(Downward  Flow) 

Eqs# 

M>1 

(Supersonic) 

M<1 
( Subonlc ) 

M71 
(Supersonic) 

H<1 
( Subonlc ) 

(8) 

dP  >0 

dP<0 

dP<0 

dP>0 

(8)&(9) 

dT  70 

dT<0 

dT<0 

dT>0 

(6) 

dOO 

dC  <0 

dC<0 

dOO 

(10) 

dV<0 

dV>0 

dV>0 

dV<0 

(It)MlO) 

df>0 

d?<0 

dP<0 

dP>0 

(12-a), 
(U)*(10) 

dM<0 

m>o 

dM>0 

dM<0 

(13).<fc) 
&(10) 

dTo<0 

dTo<0 

dTo>0 

dT,>0 

(Ik),  (k) 

8r(10) 

dPo<0 

dPo<0 

dP,70 

dPa>0 

(15), (k) 

MlO) 

dP  <0 

dP  <0 

dP>0 

dP>0 

The  Mach  Kumbsr  always  tends  toward  unity  for  upward  flow 


Developed  in  the  Appendix,  p.  73. 
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(  AZ  ;>0  ),  As  the  Mach  Number  approaches  unity,  dZ  approaches 
zero.  Hence,  there  Is  a  maximum  length  of  the  duct  for  a  given 
Initial  Mach  Number.  After  the  Mach  Number  reaches  unity,  a 
further  upward  increase  In  the  length  of  the  duct  results  in  a 
reduction  in  the  flow  rate,  i.e.  the  flow  is  choked.   If  the 
flow  is  downward  (  aZ  <0  )  the  Mach  Number  decreases  for  sub- 
sonic flow  and  increases  for  supersonic  flow.  There  is  no  chok- 
ing for  this  case. 

Analytic  Equations  for  the  Properties 

Prom  the  basic  physical  equations,  the  following  properties 
are  derived  as  function  of  Mach  Number. 


AZ*        1 

c*2  s   "k^T 


(M)   1+k    -  1 


+  -L.  (M)~k+T*  -  l 

2   k 


(16)2 


For  the  definition  of  Mach  Number,  (see  Eq,  92  Appendix) 

-        2(l-k) 
-I-  -   (M)  I+k  (17) 


With  Eq.  (17)  the  other  equations  are  derived: 

From  the  isentroplc  T-  3  relation  and  equation  of  continuity 


See  the  section  on  p.  l|l. 
Developed  in  the  Appendix,  p.  75, 


2 

-y»  -   (—7-)   ■   ("tf") 

Prom  the  isentroplc  P-T  relation 

k 

2k 

-IT  "  Nfa    L  -  <  +  > 

1+k 

11 


(19) 


(20) 


From  the  sonic  apeed  -T  relation 

1         1-k 

-4  -  <-£>  2    -  <*"  <21> 


And  from  the  definition  of  stagnation  temperature 
2Urk) 

».  .  («)  1+k      f  2     *  -M 


(22) 


Substituting  Eq.  (20)  into  the  definition  of  stagnation  pressure 

2k  k 

ji  -  (  -k- »  "^  <  -sir  -  -^-  «* »  »» 

Substituting  Eq,  (20)  into  the  definition  of  impulse  function 

2k 
I      1   _i_  %  1+k   /  1  +  kH2  »  /51.  x 

—  .  (  —  )     (  — j-5 —  )  (2k) 

The  weight  of  gas  from  any  section  up  to  the  critical  level  is 
obtained  from  the  energy  equation 


1 

Developed  In  the  Appendix,  P.  75. 
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A2        m  m  1 


p»     l-k 


?»     2 
(  •* —  I 


-  1 


Hence, 


W     .      J     ,A-|-dZ 


0 


Ar 


L    '      "  •  ii  unii    ii  (  I'      m       )  "> 


+  < 


11 ) 
1  J 


(25) 


Numerical  Example  of  Isentroplc.    Supersonic.   Upward  Flow  of  Air 

Ml  "  2»  Tl  "   ^00  R»   Pl  •  ?0°  psU»  A  "  X  f*2» 
Hence, 

Sx  •  0.5k  lbn/ft3,  Vx  -   2192. 1U  fps 

Cx  •  1096.0?  fp»,  Tcl  -  900  R 

Px  •  950U0  lbf,  Pel«  783  paia 

Using  the  energy  equation2  and  letting  Mg  *  1.75 

1/3 1       m  ( U/2.U 


AZ 


*  CUS  I1'1' 


0.2652 


34286 


+  2 


1-0.875 


AZ  » 

10t6lj.O  ft 

T2 

-  I'^hi 

2(l-k) 
1+k 

1 

-      (    1.142 

u 

ti 
•  1HJ1 

2.5 

m 

(  0.95647  ) 

2.5 


i1/3    -  1.01*55 


0.891*71 


Developed  in  the  Appendix,  p.  75. 
Developed  in  the  Appendix,  p.  7I4.. 
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3.5 


Ji    «     <l.ofc55)w    -     1.161 
pl 


W     -     p2  •  pl  ■  95,014-0  -  88,967  *     6073  lbf 

Co  .     T-    o,5 

-£•(•#»  J  -     1.0225 

Cl  Tl 

-T1    •     (  "IT  }   (     ^if*2     *     "     1*0^5  ^.69583)  -  0.9366 


-l£    -     1.1665  (0.69583)3'5    -     1.1685  (0.681)     -     0.7962 
p*l 


P2  P2 

-     <  -JT  )   (0.80111*)     -     0.9361 


This  process  was  repeated  for  other  values  of  Mg,  and  the 
results  are  plotted  on  Pig,  Ij.,  Also  plotted  on  Fig,  5  are  the 
results  for  lsentropic,  subsonic,  upward  flow  In  which  the  lni« 
tlal  conditions  are 

Hj^  •  0,5,  Px  -  100  psia,  Tx  «  500  R 


14 


o.<v- 


0.5 


Pig.   4.   Supersonic,   upward  ft-ovf1  and  subsonic   downwardlflow 
in  reversible   adiabatic  process. 


NOMENCLATURE 

A    flow  area,  ft 

a    constant 

B  notation  for  indicating  the  location  where  the  pressure  is 

the  back  pressure 

b  constant 

C  sound  speed,  fps 

C  sound  speed  at  Mach  Number  unity,  fps 

C»  constant 

c  constant 

Cp  specific  heat  at  constant  pressure,  Btu/lbm-R 

D  diameter  of  the  pipe,  ft 

F  impulse  function,  lbf 

f  friction  coefficient:  2^%/fV2 

fM  mean  friction  coefficient 

0  mass  velocity:  ^V,  Ibm/ft  -sec 

g  acceleration  of  gravity,  ft/sec2 

gg  constant  of  proportionality  in  Newton* s  second  law, 

32. 17k  lbm-ft/lbf-sec2 

J  conversion  factor,  778  ft-lbf/Btu 

k  specific  heat  ratio:   Cp/Cv 

L  length,  ft 

M  Mach  Number 

m  constant 

n  constant 


15 


o.o<j- 


0.2.O 


O.oS       o.  U        o.i6 

Fig.  5.  Subsonic,  upward  flow  and  supersonic,  downward  fl 
in  reversible  adiabatic  process. 


o.lS 


ow 
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REVERSIBLE  DIAHATIC  PLOW 

Physical  Equations 

Energy  Equation! 

CpdT     +    -Z£L    +   (  __*.    «     q  )   dZ     -     0 

-Jf-    *    J#~*      -^~    ♦  <-£-*>(  k-1  >  A-if*)  .  o 

(26) 

where  Q  is  the  constant  heat  flux  in  Btu/lbtn-ft. 

Momentum  Equation: 

jl,  +  4»«  ♦  j™l  .  0 

y     So       ge 


■¥-^  »•»  *  +  -£ 


Perfect  das  Equationt 

dP       d?       dT 


Definition  of  Mach  Number! 

dM2   m        dV2   _   dT 

]ff        y2 

Equation  of  Continuity: 
jl       1   dV2 

—  *  ^~v^ 


(27)' 


Developed  in  the  Appendix,  p.  76, 
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Analytic  Treatment 

Combining  Eqs.   (1),   (2),    (i+),   (26)   and  (27),    the  differen- 
tial relation  of  Maoh  Number  and  velocity  is  obtained.     Separat- 
ing the  variables  and  carrying  out  the  integration,   givea 


V2  -  C< 


2(   g/go  ♦  J3k  -  JO.  ) 

'  (a/so  **<ik  -W4kg/») 

JSL4-JQk.JQ+-£*  +  JQkdc-DM2" 

So  8c 


K2 


When  M  •  1,  V  •  V* 
Hence, 


V 


M^g-  +  JQk2  -  JQ  ♦  g-k) 


E  +  JQk  -  JQ  ♦  |-k  •*•  jfldnt)^ 

^  >  4 


(g/gc^J^k-JQ) 
( r/gc+ «*■&- wWSc  ? 


For  simplicity,  let 


X  -  g/Sa  •  5^  "J^  +  k8/6o 
2  •  S/Sc  ♦  JQk  -  JQ  +  kg/ge 
,  «  JQk(  k  -  1  ) 


which  results  in 
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(28Y 


Using  the  same  methods  the  following  are  obtained: 


T 
T* 


(       <(1  m2 I 


•Ci    M2 


2* 


«2  +«3  M*   * 


«,    M4 


-a  ( ir1^ 


T» 
T* 


3 

1 

"2" 


M* 


1  ♦  ^M2 


Ti 


1+k 


f      2T^      .        k-1      >    f       "l  E 


z% 


(29) 


(30) 


(3D 


(32) 


9. 
P* 


~ktr  ♦-ri 


k 

k-i_  a  >  ^^     i  .    *iM        ,  *u 

^  }  T*    «7*fJ  M*     )  (33) 


c#         "TT  v  <o  +  «c    M2        ' 


(A) 


AS 
c»2gc  "     W5^ 


Tltfl   -.  fjal-,     <  JL»  Mak  )<      "l  ** l*% 

1      (35) 


Developed  in  the  Appendix,  p.  77. 
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F*     J^d+k)     «2  ♦  ^M2 


?^ 


3  -  S1 


F  -  F* 


In 


-j|  I  ( * — -r) 


«2  + 


«t'2~ 


p7k 


( 


^«2     <¥>«* 


«c2  +  «j^ 


) 


(36) 


(37) 


(38)' 


To  determine  the  Maoh  Number  at  which  choking  occurs  (max- 
imum value  of  entropy),  Equation  (37)  ia  differentiated. 


dS 


(38-a) 


Let  dS/dM  •  0, 


M 


— £(  «jk  +  «}.  «  1  ) 


1/2 


«3 


Substituting «c2f  «c*  and*.  , 


1,2  -    X(VJf)    -  i 


This  is  the  Msch  Number  at  which  the  entropy  is  a  maximum.  Dif- 
ferentiating Eq,  (35)  with  respect  to  Mach  Number  and  setting  it 
equal  to  zero,  gives 


Developed  in  the  Appendix,  p.  79. 
■Developed  in  the  Appendix,  p,  80. 
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2  "  ~^T         1/2 
K  -  (  _- !fc ) 

*2% 
Subatitution  of  c^,  «u  and  c^,  gives 

M  -  1 
Hence, 

-$«.  -  (a/go.k)(k-i)  c ->$*-- -*p- >  -  o 

These  equation  show  that  the  length  of  the  duct  Is  at  its 
critical  maximum  length  when  the  final  *Sach  Number  reaches  uni- 
ty} at  the  same  time,  the  entropy  is  at  its  maximum  value. 

Consequently,  there  ia  a  limiting  length  of  the  vertical 
duct  for  each  initial  Mach  Number  for  the  cases  in  which  there 
are J   (1)  heat  input  and  upward  flow,  and  (2)  heat  input  and 
downward  flow,  which  is  also  the  case  for  horizontal  flow. 

A  further  increase  in  the  length  of  the  pipe  beyond  the 
critical  length  for  subsonic  flow,  produces  a  choked  condition 
in  which  the  initial  Mach  Number  is  reduced  in  value.  In  the 
case  cf  supersonic  flow,  Increasing  the  length  of  pipe  beyond 
the  critical  length  produces  a  normal  shock  within  the  pipe  or 
in  the  nozzle  whieh  feeds  the  pipe.  This  normal-shock  condition 
will  be  analyzed  later. 

If  heat  is  rejected  by  the  gas  in  its  vertical  flow,  the 
Mach  Number  is  decreased  when  the  flow  is  subsonic,  and  the  Mach 


Number  is  increased  if  the  flow  is  supersonic. 

It  is  also  interesting  to  note  that  the  maximum  entropy 
change  for  heating  with  downward  flow  is  larger  than  the  maximum 
entropy  change  for  heating  with  upward  flow,  provided  that  the 
initial  conditions  are  identical  in  both  cases. 

In  addition,  the  Mach  Number  at  the  maximum  static  tempera- 
ture is  considered.  Differentiating  Eq,  (30)  and  setting  it 
equal  to  zero,  the  result  is  obtained  as  follows. 

<?(   2-Ck  -  1  ) 

! — 

and 


«a  .  a  2> 


■  •  j  *&.**  •  y x  mM&        (3„: 


The  Mach  Number  at  maximum  temperature,  for  subsonic  flow, 
increases  toward  1/Vk  as  the  heat  flow  per  unit  mass  of  gas  (Q) 
per  unit  length  of  pipe  increases.   But,  in  the  case  where  the 
heat  input  is  small,  so  that  Q,  is  equal  to  g/(g0J),  the  highest 
value  of  static  temperature  occurs  At  the  initial  state,  and  de- 
creases thereafter.   In  the  Raylelgh-LIne  process  for  horizontal 
flow,  the  maximum  static  temperature  occurs  when  M  »  l/fc.     It 
is  found  that  the  value  of  the  tfetoh  Number  at  the  maximum  static 
temperature  is  not  a  constant.   It  depends  on  the  amount  of  heat 
influx  or  efflux  per  unit  length  of  pipe. 


developed  in  the  Appendix,  p.  8I4.. 
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Numerical  Example  of  Reversible  Jiabatlcr  Supersonic.  Upward 
Flow  of  Air 


The  initial  conditions  are 

•L 'm   2,  IW  •  100  peia,  ?x  •  $00  R 

Q  «  100  3tu/lb»-ft 
From  Eq,  (97)  (««•  Appendix  p.  77) 

HH  Ai   >^  •  (  MM  M2  )  «£»  So 


<97») 


From  Eq.  (98)  (see  Appendix  p.  78) 


(  140,2  M2  )&L  .  iC-  -  2.712  M2  3§-  g„  -  0 
V2    M2  V2   c 


(98«) 


Combination  of  Eqs,  (97*)  and  (98'),  and  elimination  of  £5»  g 

yd      m 

gives 

(  2.756  +  2.178U  M2   )  -^    «     It.1120  -^4 

M2 


dV2 
V2 


"    O.WSj+ft.gaW  W 


* 


and 


m(  -J2-  >2  -  m  f(  4s  >2<  Q.67Q?  1  o,?2?7  *?  > 

vl  I       **!  0.6702  ♦  0.5297  M2    J 


1.1*920 


Lot  K2  »   1.9 


0.9810 
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~T"  1     (  l£  ^  "7^  )2    *     (  1-lo8o>(  0.9810)2  •  1.0663 

-12  •     -_1    •  1.019^ 
?1  v2 

ySl  -     1.0663(  1  *i^2         >     -     1.0663(  0.95667  )   -  1.02 

P  9  T 

-A  «     <-_£)(«-§).     1.019M   1.0663  )  -  1.0869 
*1  yl  ll 

P        P 

r»  •  (  «JI  )(  0.656ij.  >  •  0.9306 


-^  -  ln(  1 


.0663  )  -  (  0.2857  )  ln(  1,0868  )  -  O.Oipia 


AZ  a 

• 

ft                          V      2 
(  -j2-  )   ♦  0.8(  -^2  )   -  1. 

k,        1                         1 

8009 
> 

(  jtjg  ) 

•    0.01255 

c2 

1.0326 

P2     . 

^fUjgg*  (  6.05U  )     -     0.99689 

1       . 

tfoi 

*%$^  -  ^  - 

0.003108 

This  procedure  was  repeated  for  other  values  of  Mp,  and  the 
results  are  plotted  on  Pig.  7.  Also  plotted  on  Pig.  8  are  the 
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0.04-  O.OQ 

(A2/c»jJt 

Pig.  7.  Reversible  heating,  supersonic  ^upward  flow  and 
reversible  cooling,  subsonic  downward  flow. 


0.12. 
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Fig.  8.  Supersonic  downward  flow  in /reversible  heating 
cees  &  subsonic  upward  flow  in  reversible  cooling 


-0.22 


pro 
process. 
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Pig.  9.  Reversible  heating,  subsonic  upward  flow  and 
reversible  cooling,  supersonic  downward  flow. 


~ — ' — 
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Pig.  10.  Reversible  heating,  subsonic,  downward  flow  and 
reversible  cooling,  supersonic  upward  flow. 
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results  for  heating  with  supersonic,  downward  flow.  Pig.  9 
shows  the  results  for  heating  with  subsonic,  upward  flow,  and 
Pig.  10  shows  the  results  for  heating  with  subsonic,  downward 
flow. 

IRREVERSIBLE  ADIABATIC  FLOW 

Physical  Equations 

Energy  Equation: 

-^♦¥  "^  ♦¥*!?*  -  •  m 

Momentum  Equation: 

J»  +  |_  $Z  +  **:  ig  +  kM*     fdZ     -      0  (W)1 

Perfect  Gas  Equation: 


dP  .   d?    dT 
-T    — +  T 

Definition  of  Mach  Number: 

dM2  m     dV2  m   dT 

""m2"   "v2"  "*r" 

Equation  of  Continuity: 

dS  .    1  dV2 

Developed  in  the  Appendix,  p.  85. 
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Analytic  Treatment 


Combining  the  above  equations,  eliminating  some  differen- 
tial variables,  and  rearranging  gives 


L 


eJ^JgL|«l<4l#*   (^)k(f  l^}«j 


k(  .2£  )M2  ♦  &-  .  L. 


7 


Integrating  from  the  critical  section  to  a  specified  section 

yields 

k2     2fm 
(  -g^-  )(  -g-  ) 

l/i^ 
V   I     K»<>**(fo-[.*lftM*l)jffi,) 
ln(^>  "  to  k,1c.l)(^)M2+ic(^)+[k2(i^,.(Ic2.1)|-^)1/2 

kTT 


s 

M^ 

>» 

/k+l\ 
Ut0 

S>n' 

k( 

Si 

"IT 

)   + 

k( 

¥ 

)( 

2*n  | 
TT  ' 

(►, 

•  In 


*fr4A$4fe#-  2  f*>1/2 


So 


Ttr 


A^*[i*v£/m-cAi>*-4#) 


175" 


] 


y 


lib 


I  k  So  -mz +  k( 


k  ♦  1   ,   2fm 

2       )  -5- 


Where  V  «  V*  at  M  •  1 

T  *  Mean  temperature  between  the  sections. 


f  -4- 


f  dZ 


The  details  of  derivations  are  described  in  the  Appendix  (p.  86) 
For  simplicity,  let 

Pi  "     [kVf/D  -    ^^XgSf^RT^)  j  ^ 
«      .       k(2fw/P)  -  P  1 

p2         MSTTHTfc 

^Cg/g^RT,,,) 


"TT~ 


*&-    - 


k2fn/D 

^(S/'SbSW.)   4  k(k|i)(2fB/D) 

kt^ijtafn/D) 


7  ^(s/gc2^)   +  ki^iZf^/D) 


&'J°h) 


32 


The  above  equation  simplifies  to 

v*      |Tr  *  p^  J    I  p5+  p6m2+  p7m^  J 


IHE+IT 


(U2)' 


Also,  the  other  relations  are 


JL 


P3  M2  ♦  Pli 

M2  +     p2 


M^  M2 


2fiblT 


(1*3) 


^(-^)2 


"i2 


2       «       n2P8 


M* 


£ 


^♦ftTJ       I ??  h^  +  ^  , 


(UU) 


Pe, 


i    r*2  +  P2  T°r       M1* 

a2     [  ?3-2+  fy  J      1 *$♦  P6A  V^ 


alfelT 


(1*5) 


To 

1? 


T* 


U^)    M2 


1  +  k 


2pg                       ,  k+1 

r  J*~*P2\  ,  E \  t    2**2     fc  k-1   x      n.M 


Developed  in  the  Appendix,  p.  36. 
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<4r><— - tFs > 


4  2Tk+D 

'  *E+T  *  k+1  " 


Ami*.        A      ***r.^   ■*-** 


TPpG^t  (  p54  p6m2+  PSt 


M2*  P2       ?8 


M^ 


1 

2(k+l) 


r 


F 
F» 


?*A 


1    /      5  *  P2      I    °    /    ***  v 

M2(l4k)  l^TTJ^"  „M2+  M1* 


(1+7) 
W) 

(k9) 


1f-f»1 

F»A 


^  r6' 


(  JL  )(    1+kM2   )  -   (   1+k   ) 


iilkM2!/  M^P2   i  • 


)        ( 


8T^  p3m2+  p^  >    v  T^TF^rp^sr 


1 

2(k+l) 
)             -(1+k) 

($0) 

s-s»    . 


In 


T/T** 


(^♦)<k-1,A       »8«M, 


M*+  0. 


Jt 


1 

2^ 


*     ln^{T^kfl]  (     P5+  P6M2+  P?M^   )        (51) 

*Z»  1      111     L  k-lw_M2+l2\2^8/  K^  J"*    k+1 

^  (52) 
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Numerical  Example  of  Irreversible  Adiabatlc.  Supersonic.  Upward 
Flow  of  Air 

Given  y^   «  2,   ?x   «  100  psia,  T1   *  500  R,  D  m   1  ft,  t    f  0.005 
Determine  all  the  properties  at  H^   ■  1*® 


For  any  two  sections  Eq.  (1*2)  becomes 


■  mtei 


k(k-l)(2fn/D)M|-|.(k2fn/D)«Pi     k(k»l)  (2fp/D)Mf-»-k(2fR/D)^1V 


k2  2f  m 


k(k-l)(2fB/D)M|+(k2fa/D)  +  P1     k(k-l)(2fD/D)Hf+k(2fn/D)-P1 


(k+ljDP, 


^^(g/g02RTn)+(k2fm^/D)-i-k{^5i}{2fn/D5;4       4 


^^g/gc2RTm)+(kZfwM|/D)+k{^|i)(2fH/D)M|       M^ 


k+r 


Let  Tm    •     523  *R 


px     -     0.01(  1.96-1.801/523  )1/2  -  0.01(1.956556l^)l/2 


0,013987696 


/  O.Ol8l5630li  0.050367696 


0.291932 


) 


1       lg  0.07WSJ4.029  ; 


0.20833333 


{   0.96U91itf   )(   0.97U7S068   ) 
0.91*057982 


V-,2 
Vl 


T2 


0.881t.690^0 

h£i    tip 


1S387616 

3.2i| 


1.092210 
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Tg  -     5U6     R,     Tm  -     ?0Q  |  flf6     •     523  R 

which  agree*  with  the  assumed  value  of  Tm, 
Thus  follows 

-i-S  «     -«i    .     1.06317 
91  v2 


?2     T2 


}   -     1,161205 


1A  /  y2  a 

Pi  ^ll 

CP  TP  I/2 

-5s-  -  C  «ff  1    -    1.01^51 


-p  «     ln(  -^  )  -  0.2S571U  ln(  -p2  I 
pi  1 

•     0.0862  -  O.Olj.271     •     0, 01*51*8 
Since 

5     .    UJJM    -     0.5570559 

and 

0  .  ?V    -     1183.7556 
From  momentum  equation 

2321.352  -  1*760.1^66  +  779.35512   ZK     -     0 
AZ     -     3.11*61*01*  ft 
The  weight  of  gas  is 

W»     ?„*#    -     0.5570559<0.7853982)(3.1l*61*oi*) 
-     1.37658514  lbf 
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and 
F, 


Z    -     U6t120g(l^j[l4ltlf  *  1^     .     0 .97*001*6 
|T*    -     ( 1.161205)  (0.731*33877)     -     0.85271765 


T  T 

■£*   •     {  ««jl  )  (0.915555)     -     0.999978 
°1  1 

This  procedure  can  be  used  to  determine  the  properties  of 
other  values  of  the  Mach  Number,  and  the  remits  are  plotted  on 
Pig.  12.  Also  plotted  on  Pig.  13  are  the  results  for  downward 
flow. 

NORMAL  SHOCK  FOR  UPWARD  FLOW 

Governing  Equations  for  Normal  Shock 
v2  .   "x  »  2/U-l) 


(53) 
T6X  -     T«y  <#*> 


2k       ,.2       , 


T^    "    Z—ZrTZr  (55) 


1  + 

k-1  M2 

2       * 

1  + 
1  - 

i£zi  M2 
2       7 

►  kM| 

£  •  :*3  (56) 

*x  1  ♦  kMf 
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AIL    (ff) 

Pig.   12.   Supersonic  upward  flow   and  subsoni 
In  irreversible   adiabatic  process?     SUDSOnl 


c  upward  flow 
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-u 


-IG> 


Pig.13.  Supersonic   downward  $l?w   antlf  subsonic   downward  flow 
in  irreversible   adiabatic  process. 
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Control 
Volume 


^Disconti- 
nuity 


Pig.  14.  Normal  shock  due  to 
variable  length  of  pipe  with 
a  given  initial  condition  (1) 


Fig.  15.  Normal  shock  due  to 
variable  back  pressure  and  a 
fixed  length  of  pipe  with  a 
given  initial  condition  (1). 
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Normal  Shock  Location  In  Isentroplc  Upward  Flow 

1.  Variable  length  of  the  pipe: 

M^,  P^  and  T-j,  are  given  at  an  Initial  level.  M^  la  always 
greater  than  one  (supersonic  flow),   £L  denotes  the  Increased 
length  over  the  section  where  the  Mach  Number  was  unity  when  no 
shock  was  present  (see  Pigs,  Ik  and  If?). 

For  a  given  initial  condition,  there  exists  a  corresponding 
critical  condition  In  the  pipe.  This  section  is  denoted  as  *t 
the  critical  elevation  &Z   ,  and  the  temperature  T  , 
Suppose  a  shock  occurs  within  the  pipe  at  the  height  &Xi 
The  energy  equation  for  state  (1)  to  (x)  is 


v  t*-ti  » •  \~J   *  7i-< z,  -  zi  > 


From  the  bottom  equation  of  page  7I4.  in  the  Appendix: 

f  f   M  4^         1 


•   T"   V  XT"    I  •  1 


T  I v  E7 


♦  gAK 


(57) 


The  energy  equation  for  state  (1)  to  state  («)  is 

t1 


"Vl  life 


-  1 


n?r    ,    e+i      it 


M 


giving 
AZ* 


kRT.  fe  )      1 
1    8    /    k-1 


M?    • 


*-   <lt> 


k+r  ^ 


(58) 


Since  from  Eq#  (17)  on  p,  10 

1+k 


T*  '  Tl  I  £  > 


(59) 


The  energy  equation  for  state  (x)  to  (#),  provided  no  shock 
occurs,  is 

Rearranging,  and  using  the  definitions  for  specific  heat  and 
sound  speed,  gives 

T«x     "     *£  T*  +  #•  %£<     Az*  "  Ax  >  (60) 


T4     *gc*k 
Similarly,   from  state   (y)   to  (y*) 

t-  .  ft;»if  at  ten 

Combining  Bqs.   (60)  »nd  <6l)  by  ualng  Eq.   (510,  gives 


*  *  T*-^ifc("'t 


D 


and  with  Eq,  (59),  gives 


1+k 


1  *$>     -T^ifi^'t 


(62) 


Substitution  of  Eq,  (62)  into  the  energy  equation  for  state  (y) 
to  (y*)  gives 


>•   l 


I+k 


And  using  Eq.  (53),  (57),  (58)  and  the  relation  Ay  .  Az*  + 
AL  -  AX,  the  resulting  equation  for  locating  the  shock  is 


™&c 


Ti<5^> 


2UdH 

1+k      -  A.  k-l    2    AL 


•vra 


i.( 


1-k 

4±A 


-  1 


) 


♦J 


k^x  - x 


<sf>    -  ^) 


4-     g  AL     »     0 


(63) 


Prom  the  above  equation  Mx  can  be  obtained  for  a  given  in- 
itial condition,  and  then  from  Eq,  (57),  the  location  of  the 
shook  can  be  determined. 

Numerical  example  of  normal  shock  location  due  to  variable 


kk 


length  of  pipe  for  vertically  upward  isentropic  flow  of  air. 
Given  Mx  »  2,  ?x   -  100  pala,  T1   •  500  R,  A  •  1  ft2,  find  AX  for 
a  certain  AL. 

For  convenience,  a  value  of  M  will  be  assumed,  and  the 
value  of  AL  then  will  be  determined. 
Let  Mjj  •  1.3,  Substituting  all  the  known  data  in  £q.  (63) 


(1*7007.652  -  ^) 


'  ,  0.U/2.IJ, 

~V(  1  -  I.6I883I4.      ) 


k 


,  2/2.1* 

+  JL(  l  .  0.6177285     ) 


1/3    1/3 

♦  93275(1.5381*615   -  2    ) 

k/2.k  I4/2.I4. 

♦  714-620(0.65      -  0.5     )  ♦  AL  •  0 

*L  -  UVM  •  -977-287  ft 


Substituting  M  *  1.3  into  Sq.  (57) 


37310  < j 


2  1/3 

w  -1 


+  2 


i  ,  10/6   1 
PfP  -  1 


+  AX  -  0 


AX  •  37  310  (0.6381*569)  -  23820.6269  ft 


The  results  for  different  values  of  AL  are  tabulated  as 
follows  and  plotted  on  Pig.  16 . 
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Calculation  data  for  shock  location  in  iaentropic 
flow  due  to  variable  length  of  pipe. 


Mx 

Aft 

AX 

1.0 

0 

26872.03 

1.2 

-  3^0.15 

251*56.15 

1.3 

-  977.29 

23820.83 

1.6 

-  HOteM 

15973.11 

1.8 

-  10096.63 

6683.30 

2.0 

•  16255.75 

0 

The  results  show  that  decreasing  the  height  of  the  pipe 
from  its  original  critical  height  (  AZ*)  will  cause  a  shock  to 
be  present  in  the  pipe,  provided  that  the  back  pressure  is  low 
enough  so  that  the  Mach  Number  at  the  end  of  the  decreased  pipe 
is  still  unity,  i.e.  the  exit  pressure  is  equal  to  the  critical 
pressure.  Increasing  the  height  of  the  pipe  beyond  its  original 
critical  height  will  cause  a  reduction  in  the  rate  of  mass  flow 
(choking).   It  will  be  also  found  in  section  2  below,  that  when 
the  height  of  the  pipe  Is  over  its  critical  height,  the  only 
result  is  the  change  of  mass  flow  rate  and  there  is  no  shock 
occurlng  In  the  pipe  whatever  the  back  pressure  may  be. 

2,  Normal  shock  location  for  fixed  length  of  pipe  and  variable 
back  pressure. 

In  the  case  in  which  the  pressure  is  the  only  variable  (the 
length  of  the  duct  is  fixed),  the  location  of  the  shock  will  be 
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Pig.   16.   Normal   shock  location  due   to  variable 
length  of  pipe   in  reversible   adiabatic  flow. 
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found  in  the  following  manner. 

The  energy  equation  for  state  (y)  to  state  (B)  by  using  the  M-P 

relation  Py/PB  *  (  Me/My  )2k/<1+k>  givei 


1 


• 


p„¥ 


F3 


-  1 


2 
f  PYTT 

(A    -  1 


(61*) 


The  P-M  isentropic  relation  and  Eq,  (56)  yields 


Py  - 


K 


2k 
M,  I+Ic   1  +  kM' 


1  +  kMJ 


(65) 


Substituting  Eqs.  (65)  and  (57)  into  Eq.  (61*)  with  the  relation 

AY     -    AZ  -  AX  yields 

2k 


*!n,fc'^')¥" 


2/k 


-   1 


1*^4  2^1! 

x  -4£-  &-    (    S S    )     (    5?    ) 

™r  so ■ j ;  k^i  ^ »   ^ 


1  + 


k-1 


1  +  *=! 
2 


2jk=ll 
M  1+k 

)(  5*  ) 
Ml 


kTT 


1  - 


1+k 


i  - 


(66) 


M 


Combining  Eqs.  (53)  and  (66)  gives  the  solution  of  M  . 
Prom  Eq.  (56)  AX  is  determined,  which  locates  the  position  of 
the  normal  shook. 

Numerical  example  of  shock  location  due  to  variable  back  pres- 
sure. 

Oiven  Mx  -  2,  Px  »  100  psia,  T^  -  500  R,  A  -  1  ft2,  find  Ax 
for  a  certain  P„, 

For  a  pipe  of  constant  height  Az  »  17714-2.5  ft  and  M^  » 
1.7,  Eq.  (66)  gives 

-0.28571^28  1.U2857114. 

^  (387A1168,        f  o,205U82(^8VU68) 

•  2.5  ♦  0.20511*82  +  3.21765o80  -  0.33889184 

•  2.6141132 

Hence,  two  solutions  are  obtained  from  the  above  equation: 

PB»  •  195.0  psia  (Fictitious)  &  P3  •  295.0  psia 
Also,  using  Eq.  (57) 


1/3 


k/z.h 


AX    -     37310(0.33539310)     -     12513.52  ft 


The  Mach  Number  at  the  exit  of  the  pipe,   M3  is  obtained 
from  the  isentropic  P-M  relation 

(l+k)/2k 


Vv  ■  (vpbJ 


M 


17. 


where  P  «  P^  (  1+kM*  )/(  1+kM^  ) 

"    P,    (M./M J2^1*10    <l+]df>/<l*kM*) 

x       x     x  x  y 

»  387.I4.H68  psia 

0.85711i28 
MQ  •  0.6!i055(387.ItH68/295) 

»  0,809089 

The  results  are  tabulated  as  follows  and  plotted  on  Pig. 


Calculation  data  for  the  Mach  Number  before  the  shock 
and  the  location  of  the  shock  vs.  tne  back  pressure. 


Mx 

Ax 

P3 

* 

% 

1.550 

17555.37 

353.0 

i51u0 

0,68729 

1,600 

15973.11 

33*w0 

163.3 

0.72273 

1.700 

12513.52 

295.0 

195.0 

0.80909 

1.762 

10181.88 

237,5 

237.5 

0.97867 
(  1  1  ) 

PB  Is  the  fictitious  back  pressure. 

In  this  numerical  example,  the  constant  height  of  the  pipe 
A  z  is  equal  to  17714-2.5  ft.  A  back  pressure  of  237.5  psia  is 
required  to  produce  a  shook  at  tha  section  where  M  ■  1.762, 
For  each  different  value  of  Pg,  there  is  a  corresponding  normal 
shock  existing  within  the  pipe.   If  A  z  is  reduced,  the  baok 
pressure  P^  must  be  increased  in  value  to  keep  the  normal  shock 
existing  at  the  same  location  (  i.e.  Mx  *  1,762  ),  If  A  Z  is 
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increased  over  the  length  of  177142.5  ft,  no  shock  would  exist  at 
section  x  (  Mx  *  1.762  )  for  any  value  of  PB#  In  general,  there 
is  a  limiting  naximum  height  of  the  pipe  for  a  normal  shock  ex- 
isting at  any  specified  location.  Therefore  the  "maximum 
height"  (  AZnax  )  la  introduced.  It  is  the  height  that  for  any 
further  addition  of  pipe,  the  shock  occurrence  would  be  impos- 
sible at  section  x  for  any  back  pressure.  This  result  can  be 
checked  from  the  table  on  p,  \£>  and  the  graph  on  p,  I4.6. 


AZmax(for  Mx"1#762)  •  177^2.5  ■  AZ*-  Al£  26872.03-9100 

-  17772.03  ft 


In  this  numerical  example  the  relation  of  the  "maximum 
height"  to  the  Mach  Number  before  the  shock  is  found  and  tabu- 
lated in  the  following. 


Calculation  data  for  M_  vs.  ZM._ 

x     max 


*x 

zmax 

The  "maximum  height" 

2.00 
1.80 
1.70 
1.60 
1.55 
1,00 

10782.8 
16760.0 
1931*2.1* 
21351.8 
22523.5 
26872.0 

for  Mx»l  la  Just  the 
critical  height  (  AZ*) 
of  the  pipe.  When 
AZ>26872.0,  the  flow 

Normal  Shock 

Location  in  Reversible 

is  choked. 

Diabatlc  Upward  Flow 

1.  Shock  location  for  variable  height  of  the  pipe: 
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Pig.   17.  Normal   shock  location   due  to  variable  back   pressure 
in  isentropic  upv/ard  flow. 
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The  energy  equation  for  state  (1)  to  (x)  la 

♦  «*Hi        2< 


2     x     Ml  cCg  +  e^MfT  «c  of* 


(67) 


The  energy  equation  for  state  (1)    to   (*)   Is 

Hi 


+  y^  x  ,H,  J^a&2* .  <  |-  -  »<«Mfc. 


The  energy  equation  for  state  (y)  to  (y*)  is 

(  ISL.  )  ♦  ^i(  II  )2  .  kU  f  (  ^  m   Jq  ){k.1)J^L 


^ 


Sq 


kRT 


'1 

(68) 


(69) 


y 


And  the  critical  temperature  relation  is 

2 


T*  »  T*  - 


Tk%(f7-Q,AL 


(70) 


Combining  Eqs,  (67),  (68),  (69)  and  (70)  with  the  relations 


t<2  t  c<3n2  i 


Developed  in  the  Appendix,  p.  89^ 
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T 

*7 


V         2-14.0(1,  yZ  2«|, 

1  ^  «2  ♦  •c^Mg 


and 


AY  -    AZ*  ♦  AL  -  AX 


yields 


kR 


>♦«:,  M?     2% 


y* (  ^z  x }    ■  ***?&  -  q)  al 


a>< 


o(,M 


/J i2"1*-   k-i,  «g   ,2<(^ 


k+1 


^  0(2*         Ot^My  ^        -^ 

f  M.    2-lj.«<k       «   +  «UM?     2«C^         -u  «.>+«t»M2  2o(. 

1       [**  W  *1 


♦  (£-..  jq  )(  k-1   )  AL 


(71) 


Equation  (70)  gives  the  value  of  M_.  This  value  of  M  substi- 
tuted  In  Eq,  (52)  gives  AL. 

numerical  example  of  shock  location  due  to  variable  length  of 
pipe  for  a  vertically  upward  flow  of  air. 
Olven  M^  •  2f  p  -  100  psla,  Tx  •  500  R,  Q  -  100  Btu/ft. 
Plnd  AL  for  a  certain  H^. 
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3>        *  M 


Eq.  (71)  becomes 

(  86723.961  ♦  0.22600000  AL  ) 

\     9.8688M?    1.W20171*        9.8688M?    ^<>17l* 
Sif  5.512*1*.3568m2'  *  l5.5X2+U.3568n2' 

f  M     0.981*031*8  1 .1*9201711.  M     2.981*031*8 

-  37310     A  (     22t?3?2 >  +  0.8(-S) 

I    2  5.512+1*.3568m£  2 

1.1*920171* 
{_2219J92_^) 

5.5l2+l*.3568?£ 

•  2.712  AL  -  7971*6.598 

When  Mx  •  1.2 

the  above  equation  simplifies  to 

*  1U33.6161*  -  0.0037359615    AL 

-  78515.3766  -  2.712    AL  -  7971*6.598 

Thus 

AL     -     (200.39l*)/(2.708261*l)     -     73.99352  ft 

Eq.   (67) 

1.8  •   (0.601*91327)  ( 2.700971*6)   ♦  0.8(0.21776878) (2.700971*6) 

-  (2.712AX)/37310 
AX  -  (11357. 377)/(2. 712)  «  1*187.823  ft 

The  numerical  results  are   tabulated  as  follows  and  plotted  on 
Pig.  18. 
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Calculation  data  for  shock  location  in  reversible 
diabetic  flow  due  to  variable  length  of  pipe. 


Mx 

AL 

A* 

1.0 

0 

1*615.892 

1.2 

73.99352 

1*187.823 

M 

199.13912 

321*8,331 

1.6 

281.613U5 

2137.223 

1.8 

315.13208 

1028.281* 

2.0 

315.361*35 

0 

A  normal  shock  will  exist  soujewhere  in  the  pipe  if  the  heat 
transfer  rate  is  kept  constant  and  height  of  the  pipe  is  in- 
creased provided  that  the  back  pressure  is  low  enough, 

A  special  case  for  Eq,  (71)  is  that  in  which  the  heat  flux 
Q  *  g/g0J»  AH  th«  terms  that  involve  &L   vanish  in  Eq,  (71), 
In  order  to  determine  the  location  of  the  shock  for  this  oase, 
the  momentum  equation,  equation  of  continuity,  energy  equation, 
definition  of  Mach  Number  and  the  perfect  gas  relation  are  used. 
By  the  analogy  of  the  horizontal  Fanno  Line  process  to  this 
case,  increase  of  the  elevation  over  the  critical  height  still 
causes  the  normal  shock  to  move  downward  and  finally  the  reduc- 
tions in  flow  rate  is  obtained  when  the  pressure  in  the  throat 
of  the  supersonic  nozzle,  which  introduces  gas  into  the  lower 
end  of  the  pipe,  Is  increased. 


56 


Q,ooo 


n 
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Fig.  18.  Normal  shock  location  in  reversible 
heating,  upward  flov/. 
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2.  Shock  location  du6  to  varying  the  back  pressure  with  a  fixed 
total  height  of  pipe  for  vertically  upward,  reversible  dia- 
batio  flow  of  air: 

From  the  energy  equation  for  state  (y)  to  state  (B)t 


i^  .  &     »( V^3> 


2T  "3  "B 


k-l,^,*^   ^  *2+  Sfi  x   *  j  *£.    »»%M.  ^AZ-iX 

♦  THM^      (  * 9  )       +  (ff   -  JQ)(k-l) TBW — 

■  y  "y      *2+e<3MB  7 

(72) 
Prom  the  pressure  relation: 

?B       *v  2~2"l|.   «*♦«*>£    % 

Fy  mb  'Wb 

Prom  Eq,   (56)   and  the  above  two  equations: 
1     "*     B  «2+  *3Mx 


Similarly  for  the   temperature  relation: 

,M,    *-k\  .  «,+  «,M?     2%        1  *  *,!«! 

Ty  *  Ti  (*£>  '    g      31  )  (  2_i  )  (7U> 
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Five  equations,  (67),  (72),  (73),  (7lj.)  and  (53)  relate  the 
five  unknowns,  AX,  Mj,  My,  Mg  and  Ty  for  the  given  initial  and 
independent  of  conditions,  AZ,  and  back  pressure,  ?B.  Hence 
the  location  of  the  shock  can  be  attained. 

Normal  Shock  Location  in  Irreversible  Adiabatlc  Upward  Flow 


1.  Shock  location  for  variable  height  of  the  pipei 
Prom  the  energy  equation  for  state  (1)  to  state  (♦) 

u» .  MS  %  f.  Jul  +<  x   MxS^i    )2Pe 
8  L  Mf         »A*h 


A 


l 
k+1 

) 


Ift&prf         > 


(75) 


Prom  the  energy  equation  for  state  (1)  to  state  (x) 
kRTl  8c   T     k-12    Tx   k-1  j?  vx  2 


Prom  tha  energy  equation  for  state  (y)  to  state  (y4*) 


AY*  "   kii^  I2- 


.  *+i  +(i  +  faixy  Ka ) 


A* 


2  0 


2      ' *u2    '      2   ' ' A  „2 


p3m;+  pu 


k+r 
) 


(77) 


Prom  the  critical  temperature  relation 
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7  Rk     k+1         gft 


(78) 


Prom  Eq,   (kZ) 
2 


„2  2f. 


(vS) 


2f -  o        2f  m  2f  M  2       2f  _  k+T  "TT^  1 

"k<k-l)ffi£  +k-pg  -  Pt  k(k-l)-ffi  +k~gS  ^1  >| 


-01  /» 


(79) 


T 


(5T>     T^ 


(80) 


And 


AY*     -    AZ*  +    &L-    AX 


(8X) 


Combining  Eq*.   (75),    (76),    (77),    (78),    (79),   (80)  and  (81)   and 
rearranging,   gives 

(■ 


k^TT 


H^+P2 


1 
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—     •(■ 


kR     gc 

— — -  — 

k-1    g 


♦  l^jklt.,*'1 


s+r 


pyypi,.         pT^vp^ 


V  r2 


»^-(4^)(   SP2   )       ( 1 ri 

I  n  P3*l+Pi*  P5+?6Ml+P7My  J 


UK.  6( 


•5W  <*♦¥*>♦[«£  *¥*^ 


M,    2 


k2         fm 

2f«.2     2f«  ^2       2f  Tk+ITfTir 

K  9 — ze 


2?: 


"2R 


lk(k-l)-B^*kT*  +  >1    k(k-i)^%4  +  k^.Pl/ 


¥  fc  2*1 ♦  ^ *  «¥#  4     *  k)  ***** 


^t^*^4^^  )<^) 


♦    AL 


(82) 


Bq.  (82)  combined  with  Eq,  (53) »  determines  Mx,  thua  giving 
the  location  of  the  shock. 


Numerical  example  of  normal  shock  location  due  to  variable 
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length  of  pipe. 

Qiven  J^  «  2,  ?l   »  100  paia,  T^  »  500  R,  f  «  0.005,  D  »  1  ft. 

Find  Mjj  for  a  certain  AL. 

Let  1^  «  1.2,  than  M?  -  0.7092511 

From  Eq,   (56) 


698.69(1.1279938)     -     788.11799 


From  Eq.   (Ijij.) 

T*  •  714-9.92 


From  Eq.   (78) 


M 


7^9.92  - 


Mi. 


223.86 


(83) 


Aaauraa  T.  *  714-9. 9W  °B 


Thua  T 


m 

H 
Pa 
h 

h 

H 

H 
h 

P8 


(71*9.94+788. ll8)/2  •  769.0305 

0.013991633 

0.00114.914.107 

0.16695725 

0.831453680 

0.001214.3176 

0.83229730 

0.16645947 

0.29184108 


Subatltutlng  Into  Eq.  (77) 


AY* 


T*(   0.0026260268   ) 


(84) 
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Substituting  Eq.    (83)    into  Eq.   (81*)   and  combining  with  Eq.   (8l), 
gives 

1,9693100178  -  0.0OO0117306655    AL  •  7.157213  ♦    AL 

Hence 

AL  •  -( 5.1879030 )/( 1.00001173066)   -  -5,18781*23  ft 

From  Eq,   (83) 

T*  -  7U9.92  -   (-5.l8781*23)/(223,86)  -  71*9,91*31  °R 

This  result  checks  with  the  assumed  value,  and  hence  the 
solution  is  correct. 

All  the  numerical  results  are  tabulated  in  the  following 
and  plotted  on  Fig,  19. 

Calculation  data  for  shock  location  in  irreversible 
adiabatic  flow  due  to  variable  length  of  pipe. 


M* 

AL 

AX 

1.0 

0.000000 

23.1*32873 

1.2 

-5.18781*2 

16.275660 

Uk 

-6.01*01*66 

10.225013 

1.6 

-2.750880 

6.582281 

1.8 

1.661109 

3.11*61*01* 

2.0 

6.789773 

0.000000 

The  results  show  that  for  a  decrease  of  the  height  of  the 
pipe  from  the  critical  level  (no  shock  and  M  •  1.0  at  the  end 
of  the  pipe),  there  are  mathematically  two  possible  locations 
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where  the  shock  could  appear,  when  the  exit  pressure  is  equal 
to  the  critical  value,  P*.  The  mathematical  solution  shows  that 
(see  Pig.  19)  for  &Z  "  0  there  is  a  shock  within  the  pipe  (Mx 
*  1»73)»  this  contradicts  the  original  data  of  the  problem  that 
no  shock  exists  within  the  pipe.  The  author  recommends  a  fur- 
ther investigation  for  this  shock  phenomenon.  Moreover,  there 
is  a  limiting  decrease  of  the  height  for  a  shock  existing  within 
the  pipe  when  the  exit  pressure  is  kept  at  P*.  In  this  example, 
the  limiting  lowering  length  is  about  -6.6  ft, 

2,  Shock  location  for  varying  the  back  pressure  with  a  fixed 

total  height  of  pipe. 
The  energy  equation  for  state  (y)  to  (B): 


*  •  Sty  [m  !*$»-*»  -  ¥$  >  * 

<  y  y 


(8$) 


Velocity  relations: 


2f„ 


2f. 


2f 


yB  2   f.Wl^*^  k(k-l)-^+k-32UElx 


I  k(h»l)  g*frh*pi  +  P |  k(k-i)^9*k^a  •$l) 


IT-    *-*m 


(86) 


So 


>■ 


2  2f_ 


V,  2 


k(k.i)^%2+k!£a  .  k(k.1)££aM2+k5£a  ^  u+i^  ■ 

k(k.l)^SM|+k^S  +  p  J  kik-D^^+k—S  -px 
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1.0 


3.o 


-d.o  -|.o 

AL  C  +t  > 

Fig.   19.   Normal   shock  location,  in  irreversible,    adiabati 
upward  flow.  '  x 


S.o 


*—  ,  -----  „ 
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(87) 


sc  2K1n 


Temperature  relation: 

(  ■■  )  *   (  9  )     (  |P  )  (88) 


Pressure  relation: 


(89) 


P  V  T 

(  «L  )  •  t  «M(  fi  )  (90) 

*1  vx        Tl 


And 


vy  *  vx  «  j£fe  ' 

Twelve  equations,  (53),  ($5),  (56),  (76),  (80),  (85),  (86), 
(87),  (S8),  (89),  (90)  and  (91)  relate  the  thirteen  variable*, 

TB»  PB»  VB»  MB»  Ty»  py»  vy»  V  Tx»  ?x*   Vx,  f^  and  A  X.  The 
independently  chosen  back  pressure  causes  changes  in  location 
of  the  normal  shock.  Hence,  by  selecting  Pg  as  the  independent 
variable,  A  X  may  be  found  in  fcerras  of  PB  with  the  aid  of  the 
above  equation. 
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CONCLUSION 

The  analysis  has  provided  analytical  expressions  that  can 
be  used  to  evaluate  the  properties  of  tiie  flowing  gas,  the  weight 
of  the  gas  between  any  two  pipe  sections,  and  the  location  of  the 
normal  shock  within  the  pipe  for  isentropic,  reversible  diabetic 
and  irreversible  adiabatic  processes. 

All  the  numerical  results  are  represented  graphically  so 
that  general  orders  of  magnitude  and  rates  of  changes  may  be 
found. 

Cases  in  which  the  Mach  Number  always  tends  away  from  unity 
are  (a)  the  isentropic,  downward  flow  and  (b)  heat  rejection  by 
the  gas  for  both  upward  and  downward,  reversible  flows.  For  the 
irreversible,  adiabatic  process  the  Mach  Number  always  tends  to- 
ward unity  for  both  upward  and  downward  flows,  as  in  the  case  of 
horizontal  Fanno-line  flow. 

When  the  back  pressure  is  low  enough,  adding  to  the  height 
of  the  pipe  over  the  critical  length  will  not  produce  a  shock 
within  the  pipe  if  the  process  is  isentropic.  For  reversible 
diabetic  flow,  on  the  contrary,  adding  to  the  height  will  cause 
a  shook  to  appear  in  the  pipe.  For  irreversible  adiabatic  flow, 
there  are  mathematically  two  possible  locations  of  the  normal 
shock  when  lowering  the  height  of  the  pipe  from  the  critical 
level.  A  further  investigation  is  needed  to  determine  whether 
both  answers  are  physically  possible.  Also  disclosed  is  that, 
there  is  a  maximum  height  for  which  a  shock  could  possibly  exist 
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at  a  certain  level  when  the  flow  is  lsentropic  and  upward.  The 
limiting  height  is  always  equal  to  or  lower  than  its  original 
critical  height  (i.e«  A2*) . 

The  cases  concerning  a  normal  shock  in  all  downward  flow 
processes  and  for  all  reversible  cooling  flows  whether  for  up- 
ward or  downward  velocities,  have  not  been  investigated. 
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APPENDIX 

The  Derivation  of  Equations 

#  Eq.  (8) 

Substituting  Eq.  (I4.)  into  Eq.  (5) 

£_dZ  .  £il  -  £?-  -  ££  (  5£d?-  1  ) 

-  f   (  M2  -  1  ) 

(8) 

#  Eq.  (10) 

Prom  Eq.  (3) 

Dividing  by  C  T 

t  +  Tr  kgcftf  *  (k"1)  kg^ft*  "  ° 

Introducing  the  Mach  Number 

|  ♦  ^#«$r<*M^   -    0 

Substituting  Eqs.  (1)  and  (8) 

^.if+^Kzf2  +  Mf^.1}    .  0 
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Substituting  Eq,   (k) ,   gives 


«*S*£<Ax>  *£jg  +  i£ittJl* 


T  *  5S7  7" 


i   A  fe-1*  M2 


-  1   ) 


H 


1  »   (k-1)   M2 


1  dV 
t    jl 


(10) 


#  Bq.    (13) 

Substituting  Eq,  (9)  into  Eq,  (1) 


dT  „  dP       &? 

T    T "  ? 


*   JCpdT  -X-  -  iii 


4£ 
? 


Rearranging 


*£..«<£££!.!, 


Thus  giving 


dT 


(k-1) 


if 


For 


Io    *  T(   1  ♦  £=1  M2   ) 

sfL  -  Jt  +     (k-l)M2     dM2 
TT  _nT      2+(k-l)M2  "j? 


Substituting  Eq,  (2),  gives 


3£s_.4£  ,   (k-l)M2   .  dv£_  dT  , 
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{  2—^)  *  +  .Jfcii£ 


dV 


2 


a+Oc-DH2       T       2+(k-l)M2  "v5 
Substituting  Eqs,   (12)   and  (I4.) 

flu  . 2 (k-Ddt „   ttfciil?  fli 

T°        2+(k-l)M2  1         2+(k-l)M2    ^ 

Rearranging 

dT^      tflriyU   H2  )     4J  (13) 

^T      2*(k.i)>         ?  W) 

#  Bq.   (1U) 

For  stagnation  pressure 

f.  -  p  ( 1  •  I£  m2  >  * 

dP.     dp        w?/z     a;-2 

-pT"t  +  IT^1? 

Substituting  Eqs,  (1)  and  (2),  and  rearranging 

dP*    A2  ,u  dT  l    kM2    ,  dV2   dT  , 
p<»    1    T    2+(k-l)M2    V2   * 

.  £i  +  *Li      kH2   »      kM2    dV2 
f    ™  "  2+(k-l)M2    2+(k-l)M2  V2 

Eqs.  (12)  and  (k)   give 
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It'll,  ♦  i2=2£l%Ul  ^m2       df 

■     V         2+Ck-DK2     J  2+(k-l) 


2+(k-l)M2    ? 


d?f  2k(f-M2)] 


#  Kq.    (15) 

Definition  of  the  impulse  function 

P  -  PA  +  ?AV2/gc  -  PA(  1+kM2  ) 

Differentiating 

dP  m   dP  .  kdM2 

Substituting  Eq.  (1) 

dP  ,  d?    1'     kM2  dM2 

T*  T  4  T  ♦  IS?  ** 

Using  Eqs.  (2)  and  (I4.),  and  rearranging 
dF  ,  |f   dT    kM2  (   dV2   dT  , 

T     T  +  T  *  J555  (  —2  "  T  > 

*  1+kM2      1+kM2  ? 

Using  Eq.    (12) 

i&i  1 

?  l+kM*       l+kM*   '         f        1+kM 


(Ik) 


T^        ?  *    X       t*1tm2       i*1tm2   ;         <?         i*1tm2  u^> 
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#  The  energy  equation  in  terms  of  Mach  Number  In  iaentropio  flow 
The  energy  equation  for  state  (1)  to  (2): 


cp  <  w  +  -^  ♦  ^  Vzi  >  -  ° 


v2-  if 


'   (T _-T,  )    +   'S  "*    11  +  £-  AZ  »  0 


Dividing  by  klfl^: 

T*  vf>  -  V? 


x  c     1  o  1 


Prom  the  i sen tropic  relations! 


Tl   V 


m   Uk         *    (l-k)/2 


k-1     v  k-1 


Ti       "i 


M   2(l-k)/(l+k) 


and  from  the  Mach  Huraber  definition: 


M        l|/(Uk) 


The  energy  equation  yields 


2jl-k) 


-  l 


(7-)  -  1  ♦  g^f  -  < 


jm  ■  « 


ra: 


#  Eq.  (16) 

The  energy  equation  for  state  (1)  to  (»): 
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Dividing  by  CpT* 

T*'    2  kgoRT»     k 


gcRT» 


g  *  0 


Prora  lsentropic  relations. 


h 

T 


•  k-1 


1-k 
1-k   T, 


(H,) 


w- 


(92) 


(93) 


Substituting  Eq.  (93)  into  Eq.  (92) 


2<l-k) 

■  &  % i+k 

-  1  ) 

*i( 

M^*- 

1 

) 

#  Eq.   (25) 

Since 
M  «   ( 

l±k 
q»          2 

T1 

(16) 


(91|.) 
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Substitution  of  Eq,   (9l|.)    into  (16) 


AZ*     1 
^S*  tit 


<  Jr  >     -  x 


♦  | 


Q  <fr   ' 

c  V) 


-  1 


Carrying  out  the  differentiation 


P 


m  1-k    k»2  «  2      -3 


i? 


Hence, 


W*  -  A  f  &-dZ 


*2 


I. 


f   (         1-k     k-1  2     -21 

a  I  [(f)    y     -  (?•>  j  Jd? 


o»2 


l+k       1,     J     i* 


9 


(25) 


I  Eq.   (27) 

Momentum  equation: 

dP  +  ^dz  +  m  .  o 

Dividing  by  V2 

|L.  #  J8L  <£  +    A    dvf 

v2?       Sq  v2       2g0  ^2 


Introducing  the  Mach  Number 


£  JL  +   dZ  +  JLdV;. 


P  kM2 


2  "^2" 


(27) 
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Eq.    (28) 

Substituting  Eq,   (1)   Into  Eq.   (26) 


^f^^(t^HW,|i    £* 


(95) 


Substituting  Eq.   (1;)    Into  Eq.    (95) 

iy2  t.  1        o    ,»if2 


Substituting  Eq.   (27)    into  Eq.   (96),    eliminating  dP/P  and 
Rearranging,   gives 


(96) 


.  glc^  d|       1  kM2  dvf 


1  dV2       k-1  ^  dj  _   (  |_  m  JQ  )(k.1)  M2  dj^ 


2     v2         2 


I  fkM2-l  -    (k-1)   M2]   ^ 


1-icM2  -   <  £-  -  JQ  )    (  k-1  )   M2 


S 


s 


dZ 


)?' 


Sfl  "  0 


which  gives 


o      r 


¥*-  "S • 


fi-jcM2  -(  §--  JQ  Xk-DM2 
So  sc 


(d|  gc)  -  0       (97) 


JS? 


Using  Eq.   (26) 


f  +  ¥„a^+(|_ 


J^  )(k-l)   M2   (  %,.    ) 
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Substitution  of  Eq.  (2)  into  Eq#  (26) 


fl*  .  41*  +  K|1m2  dY|  #  |i  #  jq)(k.D  i(2  (d|gc)  .  o   (98) 


Combination  of  Eq,  (9?)  and  Eq.  (98)  i'or  the  elimination  of 


(g-  -  JQ)(k-l){  ^-gi  )  -<f-  +JQk-JQ)(l  +  ^M2) 


SC 


dV2 

"v5 


+  (  JU  +JQk-JQ  ) 


dM2 


M' 


Rearranged  to  give 

dT2      2(  fc  ♦  TO  -  n   ) ^ 

V2   (  £-  4  JQk-  JQ+  &*  )  ♦  J3k(k-1)  M2   M2 
«o  6c 

Integration  yields 


In 


V2  •  - 


k<f|  +  JQk  -  JQ) 


§-+JQk-JQ  +  S-k 


c0 


Or 


In 


r  (g/gc)+«K&-jq+kg/g0  +  JQktk-DM2 


+  In  C» 
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Where  C  it  a  constant. 


V2  -  C' 


"(f-  +  JQk-W  j£-k)  +J^k(k-1)M2> 


(g/gcM<}k-:r<i+(g/g0>k 


M2 


When  M  •  1,  V  •  V*,  and 


•f-JQk  -  JQ 


V 


VT   (  f-  ♦  JQk2-JQ  +  f^fc)     1  80  Sc 


(£-  +jQk-JQ+  §-k)  +  JClk(k-X)M2 

Cc  *c 


That  is 


•(^M* 


<2*^3«2 


(28) 


#  Eq.  (38) 
Since 


S-S 

-"57 


In 


T/T* 


k-I 
UP/?»)  k 


(99) 


Substitution  of  Eqs.  (30)  and  Eq.  (31)  into  Eq.  (99) 

,M2       2« 


3-S 


M2   «0  +  «  -M2 


) 


In 


2  ^3' 
«1  M2 


j?   o(2  +«<3m2 


«V|  i£gi 


In 


k-1, 


m2    2«ir(TT,«u 


«a;*«/P 
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*  In 


( 


M2/*    «  2+  <  .M2 


jM2     (  k  **% 


(38) 


#  Mach  Number  at  the  Maximum  Value  of  Entropy,  p.  19 
Differentiating  Eq.  (36),  gives 


^M2  k%       2  1  1     k+l_        *1M*  k"V  X 


dS  _      ^2T  Zl 

c  

p 


^(- 


"i  »2    ,<¥>% 


.27^? 


(  «_♦  «<J02«<„M  -^M^*     M)     ^ 
2        3  11  3 


dM 


"  l(-|-)dM 


k+1     .. 


f    (  o(+«(M2)2«C  M-«(,M2(2«<   M)     1 
(  «2  +«(3M2)2 


dM 


-  ft  j  >«  ♦  w^HH  x  -  7^-55 '  * 


-  -  i<  |  >«  ♦  *<¥>  I «  TT71?  '  " 


.XjK2 
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Let  dS/dM  ■  0,  and 


2(1).  ilk  (  kfcl  )  ( 
kl  M  ;    M   v  k  '  v 


«C2  *«J*« 


«2  ^M2  .  x^  «2<  k+l  ) 


M2,  «2  ajL*Lii«i 


(100) 


Which  is  the  Mach  Number  at  the  maximum  entropy.   Substitution 


of  the  expressions  for  <*  t    *  and  fl  ,  into  Eq.  (100),  gives 


M2-52. 


■»■  JQk  -  JQ  +  £-k 


Sr 


r (§-  +  JQk  -  JQ)  k 


JQk(k-l) 


I 


(iL.  *.  JQk-JQ+  &-k) 
8c  «c 


♦  JQk  -  JQ  -  £-  -  JQk  ♦  JQ  -  S-k  1 


^c 


K. 


£_  *  JQk  -  JQ  ♦  &-k 


I 


!■ 


f-k  +  JQk2-  JQk  ♦  |-  ♦  JQk  -  JQ  -  §- 

m% ©o tag 


.  JQk  ♦  |  -  g* 


JQk(k-l) 


JQkfc  -  J 

35k(k-i 


-  JQk 


Differentiating  Eq.  (35)  with  respect  to  Mach  Number  yields 


<  ^2  80  } 


<§--  JQ)(k-l) 


qfM'jr« 


k-1 


)( 


•CjM2 
«2+  ^M2 


«M 
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8c     dZ 
,>*2  3R 


-1 


§*  -  JQ)(k~l) 
«e 


\   — »  +  — *- }2<<|     (      ■  ■  tH 

2        3 


2V1 


; j>  9       J +( ~z)        (-  -4)  »0 

(«2+«jr)2  *&*r 


Rearranging  and  ■implifylng,   gives 


MU-1  (j  «2+  ^3^)2  ^M-  ^^2  «3}T 


(  «2  ♦  «3H2)T 


2     (         *1m2 
M2         «  -  ♦  *  -M2 


t«I 


*2  +  V**' 


MZ  "  (  «,  4-^.M2)  U 


(    2JM 
2M2 


^2,^  2*2*^  +2*3^  -2*^  x 


)(  — = *a *—  ) 

VY  4 


(   2+kM2-M2)  ^M2 

Jfc 


1 


M2T^+^M2T"     "  % 


The  resulting  relation  is 
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2-1 

M2.        mt2M 

(101) 

* —  ♦  1  -  k 

«2  «k 

Since 

&-  +  JQk  -  JQ  +  £-k 
^           So                              Sc 

*  lj.             £-  ♦  JQk  -  JQ 

«3                                           JQk(k-l) 

«2    1       (&-  +JQk-JQ+  £-k)   g/gc  +  JQk  -  JQ 

«o                     e0       g/^  +  JQk  -JQ  +  kg/^ 

J^k(   k-1  ) 

(§-  +  JQk  -  JQ) 
8c 

g/g0  ♦  J^k  -  JQ  +(g/g0)k 
M2  .         "  6/ Sc  +  JQk  -  JQ 
JQk (k-1) 

s/^  r jS  rn  * i  - k 

(g/gp)   +  JQWQ-(g/go)k 

(g/gc)   ♦  J^  jj  JQ 

JQk2  -JQk  +(g/g;)»  JQk  -JQ-(gr/go)k-JQk2  +  JQk 
(g^gc)   +  j^k  -  n 

(g/gc)  ♦  JQk  -  JQ  -kg/gc    s 
"  (g/go)  ♦  JQk  -  JQ  -kg/go  * 


Furthermore,  when  M  «  1 


8k 


AZ»        .  1 

C^80       <JL.-j<i)(ic.i)  I 
Sc 


**{  ££)< 


1   2«, 


o<2  "  °C3 


(102) 


Since 


«1  (g/g0  ♦  JQk2  -  W  »  kg/gp) 

o(2  +oc3  "   g/so  +J*}k  -jq  ♦kg/g,,  +JQk*  -  JQk 


-  1 


Eq.  (102)  yields 


0*2°°    (&-  -  l)(k-l)   2^    •* 


#  The  Mach  Number  at  the  maximum  static  temperature  for  a  re- 
versible diabatic,  vertical  flow. 
From  Sq#  (30),  by  logarithmic  differentiation 

*1        3(2  ^3 


dM' 


•-7 


~    ,    2dM   2c<3MdM 


2  ( 


2*^  -    1      2<X  ^c< 


*  \     )   dM 


o(  24>  o(  ^M* 


Let  dT  -  0 

(2«<k  -  D(  «2  ♦«3M2)  -  2o(  «.  M2 

**  •  «  2(  2«  ^  -  1  )/« 
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That  is 


„2  .  (g/gp)  ♦*£  ■  J^  -(g/gp)*   (g/gp)+  JQk  -  J^  -(g/gp)k 
Jvik(k-l)         (g/g,,)*  J^ik  -  JQ  -(g/gc)k 

(g/gpMl-k)  ♦  J<4(  k-1  ) 


JQk(k-l) 

JQ  -  g/gp  ^z 

(    JQk    ' 

#  Eq.  tkl) 

Momentum  equation; 

AdP  +  5£-  Adz  ♦  ?AV„dV  ♦  'OfDdz  *  0  (103) 

,c         gc 

Dividing  by  A  and  introducing  fwm  ^V2^^,  Eq.  (103) 
yields 


dP+?fi-dZ4lVdV+il!^dZ.   0 

gc  e0       2^    TT4 

Dividing  by  TV2,  yields 

ii,  +  S-  dZ    dV2     1  |t  dZ  -  0 

Upon  noting  that  ?V2/gc  -  kPH2 

f  +  k*2g  *§  +  *g£  flj  *  kM2  fdZ  -  0 

M  + 1  ^  t  k^  dv2  i  d1  if  lz .  Q 


(111) 


86 


#  Eq.    (kZ) 

Combining  Eq.   (1)   and  Eq.    (Ij.)   yields 

dP  1    iV2       dT 


(10k) 


Substitution  of  Eq.  ( IOI4.)  Into  ■*«  (kl) 

^*z^  +  ^^  +  g^i§  +  kM2TfdZ,B0 


(105) 


Combination  of  Eq,  (105)  and  Eq.  (kO)  for  eliminating  dT/T 
gives 


(106) 


Also,  substitution  of  Eq.  (2)  into  Eq,  (kO)  yields 

dv2      d:;2      k-1^2  dV2      k-i  «    dZ      n 
v2  *  "J?  +  "T"1^  "^2  +  TT  i~  27  *  ° 


Combination  of  Eq.  (106)  and  Eq,  (107)  for  eliminating  dZ, 
gives 
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♦  i**2  ¥  ♦  h  b  fo 


Sc 


dM2 


Rearranging 


£:*W>  ♦*#«*♦  <¥>*$* 


Sc 


dV2 


*"      D    '  k  g0  RT    '      M2 


Separating  the  variables,  yield* 


dV< 


"    f/k+1%    g       1     i.i-2£«2.b.   ,k~l»    2f 


3R  M2 


Sinoe 


(nx+n)   dx     w     / 
(a+bx+cx^x       ./  a+bx+cx2 


ndx 


(a+bx+cx^x 


(101) 


mdx  n «~       x2 

+  rrr— ln- 


a+bx+cx2       2a       a+bx+ox2 


nb 
2a 


dx 


a+bx+cx2 
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(m.  Sk)  /     dx    +  n  ln    x2 

2a  /  a+bx+cx2   2a    a+bx+cx2 


nb, 


(*-  Sr> 


In  2c*+Mb2-4»c) 


1/2 


'(b^ac)1/2    2cx*b+(b2-Uc)1/2 

x2 


+  £ln 


a+bx+cx* 


Integration  of  Eq.  (108)  yields 


ln 


V2 


ln 


m-nb/2a 
^ox^(b2-liao)^2]  (b^Uo)1/2 
Scx+b+Cb2-^©)1/2, 


+ln(. 


n 
,2J 


a+bx+cx' 


And 


kd 


V 


k(k.i)  ^Sm2  *  ^£s  .  ^  itpS  .  (k2-!)!-  ^  5§i] 


1/2 


k(k.D  ^  *  5£ +  [k2  $ .  iaijl  ^  f£] 


20  I/2 


m 
2f. 


n 


k+T     '     IT 


(*  _J  .(k*.1)fe     1     *■] 


c       m 
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t  k  gc  2^  ^k(^T  TT'J   J 


8C 


Using  the  notations  p  lv  p  2,  fiy  l\**p   ^  and  P7»  glvet 

-_  .  < )    ( ) 

V  P|P*M  P5+P6^+P7^ 

#  Eq.   (70) 

Prom  the  energy  equation  for  state  (x)  to  (♦) 

Tox'  *fr*  +h  {   g8!-  *  )(AZ*  ~AX  I  <109) 


Prora  the  energy  equation  for  state  (y)    to  (y#) 


Combination  of  Eq.    (109)   and  Eq.   (110)   with  the  normal   shock 
property  TQX  -  TQy#   gives 

T7  m  T*  "  Txfrt;  (    gf7-Q  )AL  <70) 

#  An  Analysis  for  a  Vertical  Flow  with  Irreversible  Diabatic 
Plow. 
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Energy  aquation* 

CDdT  +  ~£  +  (  ^r  -  Q  )dZ  »  0 

f   ♦  *|1  «2  dv|  +  {  |_  WQ  )(k.1)K2  dj  go  .  Q        (26) 
Momentum  equation: 


Combination  of  the  equation  of  continuity  and  the  equation  of 
state  gives 

d?  .    1  dV2  .  dT  (lOlj.) 

Substitution  of  Eq.  (104)  into  Eq.  (1*1) 

dT   1  dV2  .  a  dZ   kM2  dV2  .  VM2  2f  *7  ,  ft         (till 

Combining  Eq,  (111)  and  Eq,  (26) 


•  0 


Rearranging 

ICM2.!)^  # 

[So 

rT 

+kJ^ 

¥• 

i  J<4)(k- 

(112) 
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Substitution  of  the  definition  of  Mach  Number  into  Eq.  (26) 
gives 

<^  .  dg  *  kjl  „2  dv£  %   {|_  .  JQ)(k.1)(^  Sc)  .  o 

(1+  ^M2)^-  -  ^4  ♦  (  t  -  J<i)(k-1)M2  4  gc  -  0    (113) 


Combination  of  Eq.  (112)  and  Eq#  (113)  for  eliminating  dZ  gives 
|(  M2-!  )(  f~  -  JQ)(k-l)  g||r 


dV2 

•3 


- (1  ♦  ¥**>  ft  * +kM*  t?  -{t  -^^  A) 


iv* 


/ 


t  ST  +ltM2  f  -<  §  -J«  '  *)    ■£•« 


Rearrangement  yields 


H--*»^y»-*t:i 


o 

^Hs^«aiiMits»<kV;*^e>4 


Jl-;-HM«-i)W  -fcy  •  Tjjjt 


ffr&-<f--J«><*-1>OT  +  i^M2 


And 
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dV2 

1* 


'k<  k-l  )£])^ 


7 


; 


Let 


^ 


2kf_ 


^3 


|tf»««    Xk-D^I-gl- 


ra     o       m 


k2f„ 


2     go 


ft.  -  JQ  )<k-l)-L-  -  =fl  -  JE=i  £-  4 


RT 


2    scRTm 


<L    -     k(k-l)  -S 


And  integrations 


•Sgi  J*  <f^/  2  <r3) 

<**♦  5  *3  5 1/2 


i 

2 
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3   4     5 


Uifc) 


Using  the  same  method,  the  equations  for  other  properties  for 
irreversible,  diabatic,  vertical  flow  can  be  developed. 
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A3STRACT 

A  one-dimensional  analysis  relating  to  the  variations  of 
the  gas  properties  along  a  vertical  pipe  with  constant  cross- 
sectional  area,  is  presented. 

Three  cases  are  considered!   isentropic  process,  reversible 
diabetic  process,  and  irreversible  adiabatic  process. 

Particularly  investigated  is  the  weight  of  the  gas  between 
any  two  sections,  the  locations  of  the  normal  shock  and  the  con- 
ditions for  producing  the  normal  shock  within  the  pipe. 

Numerical  examples  are  illustrated  and  their  results  are 
given. 


